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We generalize the idea of Vassiliev invariants to the spin network context, with the aim of 
using these invariants as a kinematical arena for a canonical quantization of gravity. This paper 
presents a detailed construction of these invariants (both ambient and regular isotopic) requiring 
a significant elaboration based on the use of Chern-Simons perturbation theory which extends the 
work of Kauffman, Martin and Witten to four-valent networks. We show that this space of knot 
invariants has the crucial property — from the point of view of the quantization of gravity — of 
being loop differentiable in the sense of distributions. This allows the definition of diffeomorphism 
and Hamiltonian constraints. We show that the invariants are annihilated by the difi'eomorphism 
constraint. In a companion paper we elaborate on the definition of a Hamiltonian constraint, discuss 
the constraint algebra, and show that the construction leads to a consistent theory of canonical 
quantum gravity. 

I. INTRODUCTION 
A. Preliminaries 



Since the early attempts in the 60 's to construct a canonical quantum theory of general relativity, a persistent 
problem has been the definition of well-defined quantum operatorial expressions for the constraints of the theory. In 
particular, the realization of the quantum Hamiltonian constraint (the Wheeler-DeWitt equation), given its nonlinear 
^-H structure in terms of momenta, has been particularly elusive. In terms of the usual canonical variables (metric and 
T-H ■ extrinsic curvature), the first impediment found was that the constraint is highly non-polynomial. The introduction of 
On , the Ashtekar variables Q seemed to provide a natural way to construct a polynomial expression for the Hamiltonian 
constraint. However, the formalism presented two significant complications: a) the variables in question were complex; 
^ , b) the polynomial Hamiltonian constraint was a density of weight two. The latter fact is crucial at the time of 
O^" promoting the constraint to a quantum operator. Since in a manifold without a prescribed classical metric (as 
I \ is the case in quantum gravity where the metric is a quantum operator) the only naturally defined density is of 
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weight one (the Dirac delta), most attempts to regularize |2|-P the Hamiltonian constraint produced operators that 
exhibited dependence on external, artificial metric structures introduced via the regularization procedure. Such 
remnant dependences on fiducial metric structures almost inevitably imply that the constraint algebra will not be 
properly implemented at a quantum level. Since fiducial structures do not transform appropriately under the action 
of the diffeomorphism constraint, the net effect is to make the Hamiltonian constraint not covariant. At most one 
could hope for formal consistency, when the algebra was computed ignoring regulators . 

An important step forward towards solving the two problems we mentioned above was a construction due to 
Thiemann [Q. He realized that it was indeed possible to construct a polynomial version of the single-densitized 
Hamiltonian constraint in terms of Ashtekar-like variables. To begin with, Thiemann uses the real connection variables 
introduced by Barbero § , thus bypassing the issues associated with the reality conditions of the traditional Ashtekar 
formulation. Barbero's variables consist of an SU{2) connection and a set of triads (throughout this paper we 
denote density weights with tildes). The basic idea to polynomialize the Hamiltonian constraint is the key canonical 
identity, 

G^^-""^- det, ' W 

where G is Newton's constant, and. 



V ^ d-'x. eabcEfE^E-e-J^ (2) 



is the volume element. In terms of these expressions the singly-densitized Hamiltonian constraint can be written 
as a polynomial operator. To simplify the discussion we focus on the Hamiltonian of Euclidean general relativity 
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(Thiemann Q has shown that an analogous construction is possible for the Hamiltonian of the Lorentzian theory, 
which involves an extra term) , 

which can therefore be written as, 

H ^ ^{AlVjFiJ-''^, (4) 

In j^] , Thiemann elaborated a proposal for promoting this expression to a quantum operator acting on diffeomorphism- 
invariant wavefunctions associated with cylindrical functions of the connection (see Q for a quick, physical, descrip- 
tion). The latter are functions constructed by considering gauge invariant quantities parameterized by multivalently- 
intersecting graphs embedded in three dimensions, usually referred to in this context as spin networks. On the space 
of cylindrical functions there are known measures [p^ that make states parameterized by non-diffeomorphically- 
related spin networks orthogonal. The proposed Hamiltonian constraint is strictly speaking only defined on the 
diffeomorphism-invariant dual of the space of cylindrical functions constructed using the measure we just referred to. 

Thiemann's proposal fulfilled several of the usual expected requirements for a quantized Hamiltonian constraint of 
general relativity. The operator was well defined and finite (even if one coupled the theory to matter) , could be made 
self-adjoint in a controlled fashion under the above mentioned measure, and was covariant under finite diffeomorphisms 
in a detailed sense. Because the operator was realized in a space of diffeomorphism-invariant wavefunctions, the 
expected result for the commutator of two Hamiltonian constraints is zero, and indeed this is the result that is 
obtained. In this approach there is no implementation of the infinitesimal generator of diffeomorphisms. Technically 
speaking, Thiemann's construction constitutes one of the first-ever consistent theories of quantum gravity. Consistency 
is a prerequisite for any physically relevant theory, but is not a sufficient condition. It is still to be understood if 
Thiemann's proposal embodies the correct physics of general relativity at a semiclassical level. 

A shadow of doubt was cast on the physical relevance of Thiemann's construction through an argument constructed 
by Lewandowski and Marolf jl^ (see also for further elaboration). Lewandowski and Marolf considered a different 
space of wavefunctions on which to define a Hamiltonian constraint. In particular, the new "habitat" (to use their 
terminology) contains the diffeomorphism-invariant dual of cylindrical functions we considered above, but also contains 
non-diffeomorphism invariant functions. Better yet, the same construction one followed to promote the Hamiltonian 
constraint to an operator in Thiemann's proposal is applicable in the new habitat. However, a crucial difference that 
stems from the correspondence with the classical Poisson algebra of constraints, is that the Hamiltonian one constructs 
should not commute with itself, but yield a commutator that is proportional to a diffeomorphism (in the new habitat 
the diffeomorphism constraint is a well defined operator yielding the natural geometric action). Yet, the Hamiltonian 
one obtains applying Thiemann's proposal commutes with itself, even in the habitat of non-diffeomorphism invariant 
states. In fact, it was also shown that several modifications of Thiemann's construction also produce Hamiltonians 
with vanishing commutators. In addition, if one artificially modified the Hamiltonians to yield a non-vanishing 
commutator, the end result did not appear to correspond to the quantization of the right hand side of the classical 
commutator 

Notice that these arguments do not prove that Thiemann's construction is incorrect in the diffeomorphism invariant 
context. It could still be the case that in that context the theory yields the correct semi-classical limit. In fact, 
Thiemann's construction in 2 -I- 1 dimensions seems to have the same difficulties as in 3 + 1 dimensions, and still 
contains within its space of solutions those of the more traditional quantizations of 2 + 1 gravity p3| . By carefully 
choosing an inner product and imposing other restrictions, the approach can be made to yield only the usual solutions 
of the 2-1-1 theory. However, it is disturbing that one does not seem to find a context where this construction 
reproduces the classical Poisson algebra with a non- vanishing commutator of two Hamiltonians. The commutativity 
might also have implications for the semiclassical limit of the theory ||lj] . 

In a separate set of developments, that has its roots in formal calculations performed some time ago in the language 
of loops and involving the doubly-densitized Hamiltonian |l5|-|l8| , another space of wavefunctions has been considered. 
These are the wavefunctions that arise when one considers in the space of gauge invariant functions of a connection 
the measure given by the exponential of the Chern-Simons action. The attractiveness of these wavefunctions stems 
from a large body of results in Chern-Simons theory that relates these states to well studied knot invariants (see 
[^9| for references), including the Vassiliev invariants. In addition, the Chern-Simons measure is in a formal sense a 
solution of the Hamiltonian constraint of quantum gravity with a cosmological constant [ pO|p^ . Moreover, the space 
is quite distinct from that of cylindrical functions in the sense that the Vassiliev invariants take non-trivial values 
on infinite sets of loops (or spin networks) whereas cylindrical functions take non-trivial values on finite sets of spin 
networks. Most of the literature on these invariants, however, has been focusing on the case of loops, whereas the 
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context that is of interest for quantum gravity requires the use of spin networks. Part of the problems we had to 
confront were to suitably generalize the Vassiliev invariants to the spin network context. Partial progress along these 
lines, but with a different motivation, was achieved by Kauffman and Lins |^ and Witten and Martin but we 
still need to elaborate on this to have an appropriate setting for doing quantum gravity. 

Why might it be interesting to consider these states in addressing the difhculties we discussed that appear in relation 
to Thiemann's Hamiltonian? Essentially because in both the space considered by Thiemann and in the habitat of 
Lewandowski and Marolf, gauge invariant operators involving Fab do not have a natural action. In spaces of these 
types, the usual way to represent such an operator would be to use the loop derivative, as is customarily in the 
non-diffeomorphism invariant Yang-Mills context. However, it was soon recognized that a conflict arises between the 
use of the loop derivative and the diffeomorphism invariance of the spaces in question. To put it naively, the action 
of the loop derivative "appends a loop of infinitesimal area" , and in a diffeomorphism invariant context there is no 
meaning to "infinitesimal area" in terms of functions of loops or spin nets. In the limit involved in the definition of 
the loop derivative, the numerator is either zero or finite and the denominator goes to zero. This is what one expects 
when one is considering what is tantamount to a "derivative of a step function" . The result should be a distribution. 
Finding the correct distribution that represents the operator is the main problem. In the new space we are proposing, 
it turns out that the loop derivative has a well defined distributional action, which can be derived from the definition 
of the invariants in terms of the Chern-Simons formulae. 

Another advantage of the kind of states we are considering is that generically they have support on an infinite 
number of spin networks. The loop representation dual of cylindrical functions has support on a finite number of spin 
nets only. This might have physical consequences. For instance, if one considers the computation of the volume of a 
finite region, the result will generically vanish for the loop representation non-diffeomorphism-invariant duals of the 
cylindrical functions (unless the volume in question contains at least one four or higher valent vertex of the finite set 
of spin networks characterizing the state). For the Vassiliev invariants, the result of the action of the volume operator 
of a finite region isgenerically non- vanishing. Thiemann already notices this drawback of cylindrical functions in 
the 2 + 1 context [|3| . Essentially, using Vassiliev invariants is tantamount to considering infinite superpositions of 
cylindrical states. To give an intuitive (and imprecise) analogy, working with cylindrical functions is like working with 
Dirac deltas in one-dimensional calculus. Working with the Vassiliev invariants is more alike to working with ordinary 
smooth functions. We will see in the companion paper that these requirements seem to lead to fruitful developments 
at least in the context of 2 + 1 dimensional gravity. 

Having at hand a definition for the loop derivative in this context, one is in a much better position to achieve a 
more realistic quantum operator representing the Hamiltonian constraint. To illustrate this point, let us recall briefly 
Thiemann's construction. It starts from the classical Hamiltonian (H) which, introducing a simplicial decomposition 
in space A, can be written as, 

H = e^^''TTih^^^h;^\^^{K^(^^),V}). (5) 

where /iq.^^- is a holonomy along one of the triangles of the tetrahedra of the simplicial decomposition and Ck are edges 
of the triangles. This expression obviously represents the Hamiltonian ^ in the limit in which the tetrahedra of the 
discretization become infinitesimal, the holonomy along the closed loop becoming Fat times the tangent vectors of 
the edges of the triangle. This quantity is later promoted to a quantum operator, essentially by adding hats in all the 
quantities involved. The resulting operator acts in the space of diffeomorphism invariant cylindrical functions. An 
important point is that in this context there is no notion of shrinking the simplicial decomposition to infinitesimal size. 
In fact, this feature allows to bypass in the construction the evaluation of an ill defined limit (the one appearing in the 
loop derivative) and to simply consider a finite holonomy as representing Fab- These steps are technically correct, but 
one can roughly expect that if one is representing Fab by a finite holonomy instead of an infinitesimal one, one would 
encounter difficulties at the time of computing commutators. A finite holonomy differs from an infinitesimal one by 
an infinite number of higher order terms, that in general will have non-vanishing contributions to commutators. 



B. Strategy 

In this paper we will discuss two main topics: a) the generalization of the notion of Vassiliev invariant to the 
spin network context; this will give a concrete footing to the wavefunctions we will consider in the companion paper 
for the action of the Hamiltonian constraint; b) the definition of a loop derivative and an infinitesimal generator of 
diffeomorhpisms in this arena, which are key to defining the constraints introduced in the companion paper and to 
the computation of the constraint algebra. 

To put this work in context, we should mention that in a previous paper p3[ | we discussed partially the generalization 
of Vassiliev invariants to spin networks for the case of trivalent intersections only. This paper expands and supersedes 
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the discussion of the previous one. In particular we will discuss several points which included incorrect interpretations 
in our previous work. In it, we mostly concentrated on the definition of the expectation value of a Wilson loop 
based on a spin net in a Chern-Simons theory. We conjectured how one could obtain framing independent invariants 
assuming factorizability, and using the notion of primitive Vassiliev invariants. Further work showed that one cannot 
really assume factorizability for intersections of valence four and higher, and therefore a small set |^ of the results of 
our previous paper do not hold in general, but only for trivalent intersections. We clarify these issues in the current 
paper. In our previous work we also started analyzing the definition of a diffeomorphism constraint, and, as a warm-up 
for our current calculations, the definition of a doubly-densitized Hamiltonian constraint. However, the final result 
was dependent on the regularization and the additional fiducial structures it introduces, and therefore there was no 
chance that the operator introduced could have the correct commutator algebra. Nevertheless, the operator had 
some appealing properties, for instance it had the expected solution with a cosmological constant (the suitable 
generalization of the Kauffman bracket), and annihilated all states based on trivalent intersections p^ ]. 

In this paper we will extend the definition of Vassiliev invariants to four-valent spin networks. The four valent case 
presents considerable technical challenges, since most literature on spin network invariants has heavily concentrated on 
trivalent vertices [ p4|j2^j2l[| . In particular the whole issue of constructing framing- independent spin network invariants 
starting from Chern-Simons theory has received virtually no attention (see pq ] for first attempts). Notice that in 
principle these are the invariants of most interest for the gravitational case, where wavefunctions are supposed to 
be invariant under diffeomorphisms. The crucial need for results involving four valent intersections stems from the 
fact that the volume operator vanishes on states with trivalent intersections. We will also see, in the companion 
paper, that the Hamiltonian constraint vanishes on the framing independent Vassiliev invariants based on trivalent 
intersections (we will, however, be able to do calculations involving the Hamiltonian constraint on related types of 
wavefunctions on which the action is non- vanishing even for trivalent spin networks) . We will also show that one can 
introduce ambient and regular isotopic invariants in this context, that they are loop differentiable, and that one can 
define a suitable diffeomorphism constraint (infinitesimal generator of diffeomorphisms) on them. We will also present 
explicit Feynman-diagrammatic expressions for the Vassiliev invariants and show that they are finite and well defined 
for spin networks. The latter definition has as a by-product that one can explicitly compute the loop derivative of 
these invariants without recurring to formal manipulations involving the path integral. 

In the companion paper we will further elaborate by defining a single-densitized version of the Hamiltonian constraint 
in terms of the loop derivative we define here, we will present several non-diffeomorphism invariant "habitats" where 
the diffeomorphism constraint is non-vanishing and we will discuss in detail the constraint algebra, showing that no 
anomalies are present. We also discuss a first application of this construction in the 2-1-1 dimensional context. 

The organization of this paper is as follows. In the next section we will briefly introduce the notion of Vassiliev 
invariants in the context of ordinary loops. In section HI we will discuss in detail the generalization of Vassiliev 
invariants to (up to four-valent) spin-networks and show that one can construct both ambient an regular isotopic 
invariants. In section IV we will analyze the loop differentiability of the invariants and the definition of an infinitesimal 
generator of diffeomorphisms on them (diffeomorphism constraint). We will summarize the results and conclusions in 
section V, and also we include an appendix where the action of the loop derivative over framing independent Vassiliev 
invariants is studied. 



II. VASSILIEV INVARIANTS FOR LOOPS 

A. Review of perturbative Chern-Simons theory for ordinary loops 

We start by recalling several results we will need from ordinary loops in order to generalize them later to spin 
networks. It is well known that in the context of loops, the expectation value of a Wilson loop (trace of the holonomy) 
in a quantum Chern-Simons theory is related to Vassiliev invariants |^,^. This expectation value can be evaluated 
in various ways. One possible strategy is to evaluate it perturbatively, taking advantage of the fact that Chern- 
Simons theories are perturbatively renormalizable |p8| . In fact, it is in the perturbative context where Vassiliev 
invariants arise in the most natural way. Let us therefore consider the expectation value of a Wilson loop W0i(7) = 
Tr(Pexp^ dy°'Aa{y)), where 7 is a loop in three dimensions ^, 



^Specifically, formulas (18) and (30) do not hold. 

■^The constant k is more commonly written as k = Aixi/k, where k is usually taken as an integer for the expression to be 
invariant under large gauge transformations. 
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< W{-i) >= / DAWAh)e 



■iScsiA) ^ 



and, 



Scs{A) = J d^xTT{A AdA+^AAAAA), 



(6) 



(7) 



This construction defines a regular isotopic knot invariant ^ £'(7, k, G). This invariant is given by a power series in 
K and depends on the gauge group G of the Chern-Simons theory in question. For G = SU (2) the invariant that 
appears is (up to a normahzation factor) the Kauffman bracket knot polynomial evaluated for a particular value of 
its variable q — exp(K). One recovers the power series if one expands the polynomial in terms of k. The coefficients 
of the power series expansion are linear combinations of Vassiliev knot invariants^ 

To evaluate the path integral perturbatively, one needs to fix a gauge and introduce ghosts. It has been established 
pof that if one considers up to three-point functions, the contributions of the ghosts produce no radiative corrections. 
We will therefore ignore them here, since at the order of perturbation we will work they play no role. Explicit 
calculations of the Feynman diagrams up to higher orders including ghost contributions are available in the literature 
in the context of loops p9| . The first relevant Feynman diagram for this theory corresponds to the two point Green 
function, 



< A<^^{x)Al{y) >= n{5"^){9a.by) = n{ 



by 



)■ 



The wavy line represents the propagator, 



47r \x - 



vT 



13 ' 



(8) 



(9) 



and the dotted line represents the invariant tensor 5°'^ . We also have the vertex. 



abc 




(10) 



The Lie algebra generators Ti'^-' (in a given representation R associated with the weight of the line going from A to 
B in the spin net) are represented by the following diagram. 



(11) 



and from now on we will drop the superscript (i?) in the understanding that it is obvious that the T's are representation 
dependent. In order to evaluate the expectation value of the Wilson loop, we start by considering the expansion of 
the path-ordered exponential. 



1 + (p dx'^Aaix) + f dx" I dy'Ab{y)Aa{x) 
y 

b 



Wa{j) = Tr 

+ idx^ I dy' I dz'A,{z)Abiy)Aaix) + 
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(12) 



If one takes the expectation value of the left hand side, one obtains in the right hand side a power series in terms of 
K. To zeroth order one gets. 



^We are using the following normalization for the Lie algebra generators Aa{x) = A^^T": Tr[r"r''] = 5"''/2, and [T°',T^] = 
^jap-yj,^ with f"^'' the structure constants of the group G. 

"^Originally, Vassiliev invariants were only considered for the ambient isotopic case, a straightforward generalization can be 
done for regular isotopy (framing dependent invariants). 
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< W(7) >^°^=dimi?G, (13) 
which is the dimension of the representation of the group G we are considering. To first order in k we get, 



where the bold face hne represents the Wilson loop 7, and, 




riia (7) 



(14) 




= Tr[r,r^]r^, 



(15) 



and, 




dx"- I dy^ gab{x,y) 



1 

8^ 



dx" 



dy^Cabc 



[x 



yf 



\x - y\ 



(16) 



The first order contribution factorizes in the product of two factor^: rn , which is group and representation dependent, 
and a^^(7), which depends only on integrals of the propagator along the Wilson line. They are called respectively the 
"group" and the "geometric" factors. The group factors are given in general by the trace of a product of Lie algebra 
generators contracted with the invariant tensors of the group. They can be thought of as the insertion at points 
in the Wilson line of the Lie algebra generators, and the evaluation of the resulting net for a flat connection. This 
diagrammatic representation of the group factors is usually called a "chord diagram". The geometric factor a^^(7) 
corresponds to the self-linking number (^(7) of the loop. The factorization in terms of a group and a geometric factor 
is a general property of the perturbative expansion of < 1^(7) >. 
For second order in k we have the following contributions. 




(17) 



with, 




Tr[r„r"r;3r^] = r2i. 



TY[TaTpT"'T% 



(18) 



(19) 



(20) 



The analytic expression of the geometric factors are obtained directly from the diagrammatic representation. For 
example, 



^We use a double index notation because at a given order there usually appear more than one factor. The first index denotes 
the order in the expansion in k, the second denotes which factor within a given order. 
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I r ry 

-[i dx"- du'' dy" dw'^gac{x,y)gbd{u,w) + c.p. ], 

^J-y Jo Jo Jo 



(21) 



where c.p. means to consider terms obtained by cyclic permutations of the first term in the indices {ax){bu){cy){dw). 

The three chord diagrams included in (^ are not all independent. They are related by the so called "STU" 
relations ||2^, 



(22) 



which stems from the commutation relation of the group generators. These relations allow to eliminate all overlapping 
diagram in terms of a combination of non-overlapping diagrams. Proceeding in this way one can define a canonical 
basis for the group factors |Q. In terms of this basis, the second order contribution can be written in the following 
way, 




<M/(7)>(')=r2ii[a"(7)]^ 



r22a^^(7), 



(23) 



where 




(24) 



In the above result two things should be noticed. The first is that, via the application of the STU relations, the 
geometric factors associated with the independent group factors will include in general a contribution coming from 
the overlapping diagrams. In particular, the combination of geometric factors associated with r2i factorizes in a 
product of first order contributions. 




Moreover, applying recoupling ^l| identities ^ one can show that 

1 



(25) 



^11 



(26) 



so the first term in ( p^ ) equals the square of the first order contribution [IF (7) >(^^]^/2 < W{"f) 

The second observation is that, as r^- depend on the group and representation chosen, this implies that the geometric 
factors 0*^ (7) embody the invariance under diffeomorphisms individually. That is, each q;'^ (7) is an independent knot 
invariant. We can also reach some conclusions concerning the nature of the constructed invariants. In principle, the 
integrals involving wavy lines have to be regularized to avoid having propagators with two of their ends evaluated at 
the same point In the context of knot theory this regularization corresponds to a framing of the knot. A possible 
framing consists, for instance, in prescribing a copy of the loop infinitesimally displaced away from the original one 
along a vector field on the manifold, and evaluating the various integrals along the separate loops. This procedure 
is not needed for all invariants. For instance, if we in detail look at 0:^^(7) we notice that in the evaluation of the 
integrals involved, the two ends of each propagator can never coincide. In the first term, one of the ends is always at 
the vertex, in the second term the ordering of the integrals prevents the ends from coinciding, if they do, one of the 
accompanying integrals is evaluated along a loop of zero length and vanishes. We refer to this case as saying that there 



^For general results about the application of recoupling theory to the factorization property of group factors, see equations 
(Q and (H) below. 
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are no "collapsible" propagators. The invariants that do not require the framing procedure are true diffeomorphism 
invariants of loops, and are called in the knot theory language ambient isotopic invariants. The invariants that require 
of the framing procedure are invariants of "ribbons" (instead of loops), and are called regular isotopic invariants. 

The explicit expressions for the higher order contributions are progressively more complicated as k increases. 
However, at any order one can apply essentially the same procedure we discussed above to isolate the independent 
group factors, and the geometric factors associated with them. The expectation value can therefore be generically 
written as, 

oo di 

<^(^)>=^^r,,a'^-(7)/«\ (27) 

1=0 i=l 

where di is the number of independent invariants at order k*. The invariant associated with < W^j) > is a regular 
isotopic invariant, that is, it is framing dependent. One can extract the framing dependence in terms of an overall 
phase factor Js^ , and be left with a framing independent invariant J{"f) defined by. 

For the case of SU{2), J(7) is the Jones polynomial for a particular value of its variable q = exp(K) and expanded as 
a power series in terms of k. We will return to this issue in terms of spin networks later. 



B. Vassiliev invariants 



Let us now show that the invariants that appear at order in the perturbative evaluation of the expectation value 
of the Wilson loop are Vassiliev invariants of order i. This is a well known result in the case of ordinary loops (see 
for instance |2^]). To do this, we recall the definition of Vassiliev invariants |^^. Given a knot invariant ^"(7), one 
starts by introducing a "Vassiliev" intersection (which we will illustrate with a black square in diagrams) , defined as, 



X 





(29) 



The invariant ^"(7) is called Vassiliev invariant of order n, if and only if when evaluated with any knot with n + 1 
Vassiliev intersections, it vanishes. One can immediately see that it also vanishes for all knots containing more than 
n + 1 intersections. This is the reason the invariants are sometimes called "of finite type" . This is also related to the 
observation that the Vassiliev intersection can be viewed as the action of a differential operator. We will see later on 
that this abstract differential operation ends up remarkably being link to the loop derivative. 

To see that the invariants we constructed via perturbation theory are Vassiliev invariants, we need certain properties 
of the path integral that can be derived by variational techniques [^,^ 16 2^j3^. For the sake of brevity we will 
not repeat the whole argument here. We will just recall that one can convert an upper in an under crossing in the 
expectation value by considering the action of the loop derivative (which appends a small loop). If one does this, one 
derives the following identity. 



< W{ 



X 



) >=< W{ 




< W{ 




K < /i(7i)T„/i(72)T" > +0(/t2) 



(30) 



The meaning of this equation is as follows: any crossing divides a loop 7 into two subloops, 71 and 72, according the 
connectivity of the loop. The invariant of the right is constructed replacing the holonomy along the whole loop by 
a holonomy divided into the subloops with insertions of the Lie algebra generators. If one now considers n Vassiliev 
intersections, one can repeat the above construction for each intersection, i.e.. 



< Wi 



n times 

A A 



)>=«;«< ;i(7i)T„,/i(72)r„, • • • /i(7fc)T"i/i(7fc+i) . . . > +O(k"+0, 



(31) 



where now 71, 72, . . . are the arcs that join one crossing with the following. The generators are inserted in each crossing 
starting from the origin. Notice that each crossing is traversed twice, and that the two generators associated with a 
given crossing are contracted between themselves. We see that the left hand side of the above equation has vanishing 
coefficients up to order n, that is, the coefficients of the expansion up to order n behave like Vassiliev invariants of 
order equal or less than n. A similar behavior is observed for the case of spin networks. 



III. VASSILIEV INVARIANTS FOR SPIN NETWORKS 

A. Diagrammatic notation for the expectation value of a Wilson net 

A spin network is constructed considering a graph F (with multivalent intersections) embedded in a three 
dimensional manifold. The graph is composed by a set of edges {ej}, (j = 1, . . . ,N(r)), which are analytic oriented 
open paths that intersect only in the initial and final points. The intersection point of three (or more) edges defines 
a tri (or multivalent) vertex v. We denote N{r) the number of edges, and V{r) the number of vertices of the 
graph. Associated with each edge in the graph is a representation labelled by Ji of an element of the gauge 
group G. For instance, given a connection, one can construct these elements by considering the holonomy along the 
given edge in a given representation. The group elements are "intertwined" into gauge invariant objects through 
contraction at the vertices of the network with the invariant tensors of the group. Given a graph F and an assignment 
J = (Ji, . . . , JN{r)) of representations to each edge plus a set of intertwiners I = (Ji, . . . , Iv(r)), one can construct 
starting from a connection an object we call "Wilson net", Wa{s) = Wa(F,/, J), in the following way. 



Wa{s) 



n 



Ci...C„ 

vBi...B„ 



(32) 



We will use the compact notation s denoting the graph, intertwiners and spin weights together. In the above expression, 
all the matrix indices of the holonomies U{e) are contracted with the indices of the intertwiners 1^ in a way which 
depends on the connectivity of the net. The edges of the net are oriented according to the following conventions, 

5„ 



(33) 



and, 



J Ci...C„ 




(34) 



The invariant tensor /, 



Ci. 

vBi. 



■Bm 



can be expressed in a basis of n„ 



— 2 Clebsch-Gordan intertwinersQ. This 
construction is feasible for any group, but we will concentrate in the following sections on SU{2) and the representations 
will be labelled by integers. A more elaborate labeling is needed for other groups. Spin networks were originally 
introduced in the SU(2) context by Penrose | |3^ . 

We now wish to compute the expectation value of a Wilson net in a Chern-Simons theory. 



< Wis) >= / DAWAis 



E{s,k) 



(35) 



'^In general, a n + m-valent intersection can be thought of as the limit where the intermediate lines are of vanishing length, of 
n + m — 2 trivalent intersections connected by additional links, in a choice of basis such that the intermediate links are assigned 
a definite spin value. 
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This integral has been studied by Witten and Martin and KaufFman and Lins |21|, estabHshing the crossing rela- 
tions (skein relations) and recoupling relations for networks including up to trivalent intersections. Their calculations 
were based on the use of monodromies and the tangle group, respectively. For trivalent vertices we have also shown 
in a previous paper |25[| that one can extract the framing factor and construct ambient isotopic invariants. 

Here we would like to study if analogous results can be found in the case of four-valent intersections, and also 
perform perturbative studies of the path integral with the aim of constructing explicitly ambient invariants. We 
therefore proceed as in the case of loops, to consider order by order the perturbative evaluation of expression (|3^). 
The general form of the perturbative expansion can be easily deduced from the definition (32). We first introduce 
a convenient diagrammatic notation. The holonomy associated with a generic edge e is written in the following 
form, introducing a diagrammatic notation for the explicit expression of the path ordered exponential in terms of the 
connection, 



C/(e) 



(36) 



where, 



Be 



(37) 



and. 



J e Jo Jo 



)...5{xi -yi). 



(38) 



For the readers familiar with the extended loop formalism ( ||T^ chapter 2), the above multitangent corresponds 
to X(e)"i^i '""^" in such a notation. In the above expressions, we have adopted the usual "generalized Einstein 
convention" of extended loops, in which a summation is implied by any repeated index Oi and an integral over the 
three manifold is implied by any repeated index Xi. 

Some general results about the factorization of the geometric factors (see equation (p5|)) can be derived from the 
following property of these fields. 



E 

P(l...m,m-t-l...n) 



Op X a X 



(39) 



where P(l . . . m, m + 1 . . . n) is any permutation of the indices which preserves the ordering of the subsets (1 . . . m) 
and (m + 1 . . . n) between themselves. This equation express the "algebraic constraint" of the multitangent fields, 
first studied in j39j. Introducing ( p6| ) in (|3^), we get the following diagrammatic expression for the expectation value 
of the Wilson net. 



< 



wa{s)> = E ■■■ E 



A=flat 



m— nj\ 



^ < (^aixi ■ • ■ ^ar^lxr.-^ ) ■ ■ ■ >, 



(40) 



Notice that in the above equation one takes a sum for each edge of the net. In general, WOi=flat(s) defines the 
"chromatic evaluation" of the spin network s. This quantity is evaluated replacing in the functional integral all the 
holonomies by identity matrices, U{e)^^ — > 6^" for all e. The resulting expression can be computed using recoupling 
theory. Then, in (^ ) M^.4-flat( . ' ■ ■ ■) has the following meaning: take the chromatic evaluation of the spin network 
replacing the holonomy U{ei) by a product of Lie algebra generators Ta^ . . . T^^ in the representation if Ui ^ 0, 
or by the identity matrix if = 0. 
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At this point it is that the convenience of the notation we introduced becomes clear, we see in equation (40) that 
the expression for the expectation value of the Wilson net becomes really compact, involving a chromatic evaluation of 
a diagram times another diagram contracted with the propagators. This is quite analogous to what happened in the 
case of ordinary Wilson loops. To generate from (^) the perturbative series one has to introduce the Feynman rules 
to express the Green functions < {A'^^^_^ . . . Aal^^x^-^ ) ■ ■ ■ > in terms of wavy lines (the propagator and the vertex)^ 
and dotted lines (the invariant tensors of the group). We immediately see that a mechanism similar to that observed 
in the case of loops operates here: the contraction of the invariant tensors with W^A-flat( ' . ■ ■ ■) generates dotted 
diagrams (the group factors), whereas the contraction of the two and three point propagators with the multitangent 
fields generates wavy diagrams (the geometric factors). We schematically represent this procedure in the following 
way. 



/ 



WOl=flat 



7" (dotted diagram). 



(41) 



and, 



{gaixia„-^x„^ ■ ■ ■) — I :2 



(wavy diagram). 



(42) 



These relations define the group and geometric factors for spin networks. Notice that in the above examples we have 
chosen a specific way to connect the indices ai,ani and aiXi, amXni- In general, the application of the Feynman 
rules produce, for a given number of indices, a superposition of diagrams with different connectivities of the dotted 
(or wavy) lines. Therefore, we can schematically write (ignoring the ghosts as we discussed before). 



< Wis) >-- 



rn—0 {diagrams} 



^ ^ r[ (dotted diagrams) ] (wavy diagrams) '"^^ 



(43) 



The index m gives the order of the perturbative expansion, and for a fixed m we sum over all the diagrams that 
contributes to this order. Each term factorizes in the product of a dotted diagram (the group factor) and a wavy 
diagram (the geometric factor), both having the same connectivity (i.e., in both diagrams the connection of the dotted 
and wavy lines are identical). 

In principle, the group factors can be evaluated using recoupling theory. For SU{2) one uses the following version 
of the Fierz identity (the standard identity between the Lie algebra generators and the Clebsch-Gordan coefficients 
in the fundamental representation). 



y(J)y(A') ^ 



K J 

= (-1)2( W+^A.K- 



(44) 



where Kjk = AjA^f and Aj = ^JJ{J + Yj{2J + Y) . The Fierz identity allows to rewrite any dotted diagram in terms 
of the chromatic evaluation of a new spin network where the dotted lines are substituted by edges of spin one. In 
general, the group factors satisfy some identities which are characteristic of spin networks and which play an important 
role at the time to identify the topological invariants associated with the perturbative expansion. To introduce these 
identities, let us consider a generic dotted diagram and let us isolate from it a single dotted line connecting two edges 
e and e'. We draw this single line enclosed into a circle. We also draw one of the vertices of the selected edges, for 
example, the origin of e', which we assume to be four-valent. The dotted lines outside the circle remain fixed as we 
move the generator starting from e' all around the vertex. Then it is easy to prove that. 





+ r 





(45) 



As in the case of loops, we will ignore the ghosts. It one wishes to include them, one has to incorporate in the analysis the 
ghost's propagator and the ghost's vertex. 
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where r is any of the group factors we have considered (we omit indices since the expression is vahd at any order). For 
a given diagram, one can perform these movements with aU the dotted hnes of the graphs producing a hnear system 
of equations for the group factors (of course, not all the possible movements generate independent equations). This 
set of identities is valid for any group, and they are rooted in the gauge invariance property of the Wilson net. For 
this reason, we call these identities the Gauss' law for spin networks. 

Another set of identities, called "STU relations" arise from the generator algebra, 

^ ^ ^ + 

" Y 



The STU identities, together with Gauss' law, form a set of linear identities that can be solved order by order in 
the perturbative expansion, allowing to find all group factors in terms of a reduced set of "independent factors" . 

The zeroth order contribution of (^) is given by the chromatic evaluation of the net, which we usually denote in 
the form, 

< Wis) >(°)= WA=a.t{s) = E{s, 0) (47) 

In the following subsections we are going to study the detailed structure of the group and geometric factors for the 
first and second order in k. 




B. The first order contribution 

Considering equation (|4^) to first order in k we have, 

< Wis) >(^)= J: J: r f _JT\^ -EE (48) 

where by i and j we denote any pair of edges of the spin network (they can also be the same edge) . We usually use 
the compact notation for the first order group factors, and tpij for a wavy diagram connecting the edges and Cj. 

We notice a difference with the case of loops. In that case there exists only one group factor, , at first order. Here 
we have a number of first order group factors equals to N{T){N{T) + l)/2. These quantities are not all independent 
due to the Gauss constraints (^5|). The Gauss identities generate for a given graph a linear system of equations which 
allows to express some of the in terms of a small subset, which we call the set of "independent" group factors 
(IGF). The determination of the IGF allows to construct the topological invariants associated with the spin network. 
As in the case of loops, the linear combinations of geometric factors that multiply each group factor are (regular) 
knot invariants. 

It is clear that the IGF set is not unique (we can choose different free parameters to express the solutions of the 
Gauss identities). 

Let IGF = {zi}, i — 1 . . . o?(s), be the free parameters with respect to which the linear system of Gauss equations is 
solved. The number d{s) of parameters depend on the particular spin network being considered. Then, ( [48|) can be 
rewritten in the following way, 

D{s) 

< W{s) ^^^^('^) (49) 

j=i 

where the 2'((/9)'s represent knot invariants formed by the linear combination of geometric factors associated with each 
independent parameter. 

Let us illustrate the procedure with an example. For the net shown in figure |^, d{s) — 6 and we choose IGF = 
{ri2,ri3,ri4,r23,r24:,r34,}. Then one gets 

-?'12(</?11 + -^22 ~ ^12) - ?'34('/'33 + fiA - fSi) 
2 4 

+ E E '^^i ^ "^55 + </Jj5 +'Pj5 + ^ij )■ (50) 

i=l j=3 
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FIG. 1. A spin net with a four-valent and a trivalent vertices. 

The invariants have a simple geometric interpretation in terms of Unking numbers of loops constructed with the edges 
of the graph. It is immediate to see that each of the above linear superposition of geometric factors combine to give 
the self-linking number of a closed path 7 with support in the net, 



< W{ 




^ > 



i=l j=3 



with0, 



(51) 



(52) 
(53) 



For any other IGF, the invariants are written as a combination of the above self-linking numbers. Notice that in this 
example all the first order invariants are framing dependent. 



C. The second order contribution 



Considering the contribution to second order in k in ( |4^ ) we have. 



<M/(.)>(^)=^^^r 

i j k 



I j k I \ 





(54) 



''^X A. 



Guided by the case of ordinary loops, we would like to study if it is possible to extract a pre- factor that is dependent 
on framing, therefore isolating a framing-independent knot invariant in this expression. In the case of ordinary loops 
the exponential of the first order contribution was the prefactor that allowed to isolate the Jones polynomial. Here we 
would like to check if this is possible for spin nets up to second order in the expansion, for four valent intersections. 
We will see that, as we independently proved in our previous paper, it is possible to extract the prefactor for trivalent 
vertices, but it is not possible for four and higher valent intersections. This was already anticipated in the work of 
Labastida, who noted that it does not hold for multiloops, which are a particular case of spin networks. 

Using the Gauss and STU identities and the algebraic constraint we discussed before, and observing that. 



1 



(55) 



Overbars denote edges traversed in opposite directions, and the o operation is the usual composition law of open paths. 
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and introducing the "shifted" four-point group factors f, defined in the following way, 



it is found that, 

< W{s) >(2) 



1 



(56) 



2E{s,0) 



[< W{s) >(i)]2 



■EE 



-[ifijf +r 




E^ 



z'^ + EEE'^ 



fijVjk + cyclic permutations in i,j,k 



(57) 



i j k 

The quantities (pij are the first order geometric factors linking the edges and Cj, and. 



Pi 



+ 



(58) 




(59) 



and. 



Pijk 




(60) 



Looking at equation (|57| ) we see that we have not achieved the desired result. If we look at the first term in the right 
hand side of this equation we identify the expansion to second order of the exponential of the first order correction, 
exactly like in the case of loops. However, there remain framing dependent contributions in the other terms, therefore 
one will not be able to extract the framing dependence in an overall prefactor, as happened in the case of loops. The 
additional framing dependence arises for instance in the terms involving Lpij or pijk when the edges /c share a 
common vertex. 

A remarkable exception happens in the case of trivalent intersections. In this case the group factors r that appear in 
the above mentioned terms vanish identically. This can be seen straightforwardly applying recoupling at the vertex. 
In the case of trivalent vertices this yields only a prefactor and the two terms involved in the definition of the f, 
equation ( p6|) vanish. In the case of higher order vertices, the application of recoupling leads to a sum of terms and 
therefore the two terms in ( p6| ) are different. This fact will also appear at the root of our proof that the loop derivative 
annihilates trivalent framing independent invariants, as we will discuss later on. 



D. Topological invariants for spin networks 

The analysis of the higher order contributions can be performed using the same techniques we have developed at 
first and second order. The results can be summarized as follows. To any order in the perturbative expansion, there 
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exist a number of independent group factors for a given group G, zji^\ j, I, G), based on solving the Gauss and STU 
relations. In terms of these the perturbative expansion can be written as, 

<W>^^^Y1 ^^"^ G)^-"^ (61) 

where the la"^ are (regular or ambient) isotopic invariants of the embedding of the spin network F, and d"{s) is the 
number of independent group factors at n — th order for a given spin network s. This is a difference with the case of 
ordinary loops, where the expressions are given once and for all loops. This is related to the fact that the spin nets 
are a much more general object. 

The regular or ambient isotopy properties of the invariants can be analyzed systematically, and, to any order in k, 
the invariants are grouped in two categories: framing independent and framing dependent. The framing dependent 
part cannot be exponentiated like in the case of loops. The general result can be written in the following way, 

oo rff (s) oo 4(s) 

< W{s) >= E E ''ail I, G) lliiT) + 'fbiJ, I G) 2* (r) ,\ (62) 

i=0 a=l 1=0 6=1 

In spite to this dependence with the particular spin network considered, the invariants have the following property: 
given an embedding F, let F' e F be a subgraph obtained by removing some of the edges of F. Then, the invariants 
2'(F') are obtained from T{T) putting zero all the wavy diagrams with legs attached to any of the suppressed edges. 

The specific expression of the invariants Ifd and Ifi that appear in expression ( |62|) and their quantity d{^{s), 
depend on the choice of independent invariants picked to solve the constraints among the group factors. We will give 
a procedure that generates, for a given spin network, a basis of independent invariants that are framing-independent. 

A general invariant of n — th order is in general given by a linear combination, 

J(")(F)= ^ c„l(")(F), (63) 

a=l 

where the 2a 's have the following expression, 

2(")(F) ^ ^fa«!"^^r) + Y: Or^at'^^'in (64) 

i i 

where the D^^"^ are well defined numbers that depend on the connectivity of the spin net and the a's are products 
of multitangents contracted with the geometric portion of the Fcynman diagrams (propagators and vertices) and can 
either be framing dependent or independent. We can therefore conclude that an invariant will be framing independent 
if the choice of coefficient Cq's is such that, 

a=l 

One now solves this system of equations, which will lead to expressing some of the Ca's in terms of a set of independent 
Cd's. For each of these independent Co's one has an independent invariant. 

Let us consider a concrete example of this procedure. We will carry out the computation for the spin network of 
the figure (|l]). The invariants that appear at second order in the calculation have the generic form, 

2 4 

2(2) (D = cip(7i2) + C2p(73i) + E E ^nPiln^), (66) 

1=1 j=i 

where 7^ are the subloops of the spin network we defined in equations ( |52| , p3|) . The quantities p are closely related 
with the second coefRcient of the Alexander-Conway polynomial |2|] . These were the expressions one directly got 
when one performed the calculation in the case of loops. In the case of spin networks one gets combinations of 
these invariants evaluated in sub-loops that form the spin network. The expression is manifestly framing independent 
since the coefficients of the Alexander- Conway polynomial are framing- independent. We see that (at least up to this 
order in the perturbative expansion) the use of spin networks does not introduce new topological invariants, but just 
combinations of invariants evaluated in the various subloops that form the spin network. 
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Finally, we show that the invariants we have obtained are the natural generalization to the spin network context 
of Vassiliev invariants. To see this, consider an arbitrary net, and construct the Vassiliev intersection. Following the 
same type of argument as we did in the case of ordinary loops one immediately obtains, using variational techniques. 




= K<... h^'^\e^,)T^J^^h^'^\e2,) . . . /^(•''=He3fc)Ty'-)/i(^^)(e4fc) . . . > (67) 

where the indices in parenthesis refer to the representation in which one is considering the holonomies and the Lie 
algebra basis elements. We therefore reach the same conclusion as in the case of loops: an invariant with n Vassiliev 
crossings is of order k" meaning that the first n — 1 coefficients in the series expansion vanish for spin networks with 
n crossings. Therefore, they are Vassiliev invariants of order up to n — 1. 



IV. THE DIFFEOMORPHISM CONSTRAINT 



A. Loop differentiability 

We now have succeeded in setting a suitable "arena" where we can start discussing the constraints of quantum 
gravity. The arena will be the space of Vassiliev invariants, and in order to implement the constraints as quantum 
operators we will take advantage of the fact that these invariants are loop differentiable. This is a result we introduced 
in a previous paper p3t, we only recall the appropriate formula here. The result appears simply by applying the 
definition of loop derivative [0 to the path integral defining the expectation value of the Wilson loop. One acts with 
the loop derivative on the Wilson net in the path integral. One obtains as a result a Wilson net with an Fab inserted 
at the point where the loop derivative acted. One then re-expresses the Fab m terms of a functional derivative of 
the exponential of the Chern-Simons form with respect to the connection. This functional derivative is integrated by 
parts to act on the Wilson loop again. The result is the expectation value of a Wilson net with the insertion of two 
Lie algebra generators. 



^ab{K)E{' 



-2k ^ dabc / dy'^S^ix -y) X 
k •'^k 



(68) 



where Vj and v'^ are the start and end points of edge Cj. The loop derivative depends on a path tt^ going from the 
basepoint o to a point x, and this operator acts on the space of spin networks with a fixed basepoint o G F. The links 
involved in the formula are shown in the figure. Using (Eh one can rearrange the above expression in the following 




■J 

FIG. 2. The basic configuration that defines the loop derivative. 



way. 



o,«)==«:^(-l)2(-^^+'^^)Aj^.j,e,fc, / dy'S'ix-y)E 



A 

o 


% 1 
1 — »^ 






V ^ 


X 


'y 


, K 



(69) 
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We now replace the two parallel lines in terms of a sum of irreducible representations. The resulting expression can 
be cast using recoupling identities]^ as, 



k 



e 

J 


I 


O 1 


K 


V ' 





(70) 



For future reference, we would like to introduce a formula for the loop derivative when the path tt added is of 
infinitesimal length. This is useful in the definition of the Hamiltonian constraint. It could arise when the point x 
is infinitesimally separated along the same edge as the origin o or it could correspond to two infinitesimally adjacent 
edges (for instance incident on a vertex) . More precisely, in both the diffeomorphism and Hamiltonian constraint one 
has the loop derivative evaluated along an infinitesimal path multiplied times a regularized Dirac Delta function. In 
that case one can rearrange the above expression in terms of a chord diagram. 



lime^ ' -"^ 



o, k) = lim e^^y^ Cahc / dy''xe{o,y)E 

e^O ^ — ' /„ 




(71) 



where x is a (symmetric, Xei^^y) = XeiUjx)) regulator such that lim^^Q Xe{x,y) = 6^{x — y). We will derive in the 
appendix the rather remarkable result that the kind of derivative we considered in the previous formula annihilates 
the framing-independent Vassiliev invariants if the spin networks consider have trivalent or lower intersections. This 
will allow us, in the companion paper, to find large families of solutions to all the constraints of quantum gravity, 
albeit with vanishing volume. 

Notice that formula ( ]7C| ) implies, if we analyze it order by order in k, that the loop derivative of a Vassiliev invariant 
of order n (the coefficient of k" in the expansion) is proportional to a Vassiliev invariant of order n — 1. That is, the 
loop differentiation operation lowers the order of a Vassiliev invariant by one. 

As an instructive example of the consistency of the formalism, we can compare the result predicted for the loop 
derivative of the first order term in k by equation ( [70| ) and the direct computation of the loop derivative of expression 



Ait^«)<M^(.)>(^)=-^Ail)«)( E / f ^'''^'^f^^J^. 



(72) 



where we denoted that the loop derivative acts along the edge by a superscript. The loop derivative will only act 
on the loop dependence of the right-hand-side, i.e., on the integrals. To evaluate this action, we recall the well known 
resuh (i.e. |l|]). 



Ait^(7r:) / d2/'=F,(y)=5[„F,](a;), 



and therefore get. 



Ai?^«) < W{s) - I dz%aemdf^^ 



(73) 



(74) 



which we can rearrange as, 



AttIx — z\ 



(75) 



The recoupling relations reduce the sum over the irreducible representations to a single term with spin one. This is a 
difference with our previous result -equation (28)- of reference |^^, where we used a slightly different version of the Fierz 
identity. The one we are using here is more useful in calculations. 
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The last term gives the integral along an edge of a quantity differentiated with respect to the integration variable. 
This gives contributions only at the ends of the edges, and it gives the same contribution per each incoming edge. 
Recalling that the summation of all r^j's for each vertex vanishes, the last term vanishes when added along the whole 
spin network, vertex per vertex. The remaining term gives the final result, 

^it\<) < Wis) J2 ^^^''^"^ I dy'^'i'^ - y)- (76) 

This expression coincides with the right hand side of ( [70| ) evaluated to first order in k. The group factors are given 
in (70) by the chromatic evaluation of the invariant E, which is what one gets at zeroth order in n. 

It is worthwhile emphasizing that although expressions like ( [70| ) can be viewed as "formal manipulations" in the 
path integral that one uses to motivate certain results, yet the perturbative Feynman-diagrammatic calculations 
yield order by order results that are consistent with the formal manipulations. We will not discuss in detail similar 
calculations for invariants of higher order, but they can be checked in detail using the formulae for loop derivatives of 
multitangents, see for instance |p7| . 

An interesting observation is that if we go back to expression (|68|) we see that the loop derivative introduces in 



the spin network two Lie algebra generators contracted. This is precisely the same contribution we got in (67) for 
the Vassiliev derivative. This points out to a very appealing connection between the loop derivative and the Vassiliev 
derivative. It therefore highlights why Vassiliev invariants play a special role among knot invariants from the point of 
view of quantum gravity: they are naturally "differentiable" knot invariants, and the derivative in question is nothing 
else but the loop derivative. This is consistent with the fact that the loop derivative "lowers the order" of a Vassiliev 
invariant, as we saw before in the example of the invariant of order one being related to the zeroth order invariant. 
Freidel was the first to notice, and is currently studying the implications of this connection ||40(|. 



B. The diffeomorphism constraint and the loop derivative 



To define a diffeomorphism constraint we will formally integrate by parts the expression of the constraint in 
the connection representation acting on a state given by the loop transform. In terms of the new variables, the 
diffeomorphism constraint is given classically by C{N) = J cPxN°'{x)E\FI^^. If one wishes to promote this to a 
quantum operator one needs to choose a factor ordering and a regularization. We will choose the factor ordering in 
which the triad is at the left. In this factor ordering, the operator formally fails to generate diffeomorhpisms on the 
wavefunctions in the connection representation. This can be amended if one regularizes the operator by point-splitting 
the operator E and F and using a symmetric regularization ||l^ . We will therefore operate with such an operator on 
a state in the loop representation defined via the loop transform, 

C{N)'^{s) = j DA{C{N)-^{A))Wa{s) (77) 

where again x is a (symmetric, Xe(x, y) = XniUi ^)) regulator such that limg^o Xc(a;, y) = 5^{x — y), and the symmetric 
regularization refers to taking the symmetric average of the vector N at the points x and y. This expression is not 
gauge invariant due to the point splitting. This can be fixed by introducing small pieces of holonomies along paths tt^ 
that connect the points x and y and inserting them to get a gauge invariant expression. If one now formally integrates 
by parts the functional derivative and represents the Fab using the loop derivative, we get, 

C{N)^{s) f DA'f (A) lim j d?x j dy''^^^^^M±^^^x.(x, y)A,,(^^)W^^(s). (78) 

We therefore identify the following operator as the infinitesimal generator of diffcomorphisms in terms of spin nets, 

C{N)^{s)^Y.HI '^'"'l + xe{x, y)Aafc(7rg)v&(.). (79) 

The last expression should be understood as a shorthand notation for ([78|). This means that the limit should be 
taken before evaluating the functional integral, otherwise we would fall in the usual pitfall of not having a defined 
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action for a loop derivative on a knot invariant. This action is defined by formula (|7^) for the invariants in question, 
replacing ^'(A) by the exponential of the Chern-Simons form. The last expression has been well known for some 
time to correspond to the infinitesimal generator of diffeomorphisms on Wilson loops evaluated along smooth loops. 
In the case of spin networks, the expression displaces each edge infinitesimally. The vertices are left fixed, but are 
reconnected to the shifted edges by retraced paths. Due to the gauge invariance of the intertwiner operators, the 
functions behave as if the vertex had been shifted along with the edges. 

In our previous work in terms of loops (l^,^ we have used an (unregulated) diffcomorphism constraint that looked 
slightly different than the one we consider here, but that coincides in the limit, when acting on functions on which 
the loop derivative yields a smooth result, like for instance holonomies of smooth connections. The main difference is 
that in the regulated operator the path tt extends out of the original loop and tends to it in the limit, whereas in the 
unregulated expressions one simply took the path directly along the loop. 

The formul a (fz^ ) involves evaluating the loop derivative for an infinitesimal path. To evaluate this, it is convenient 
to go back to (|6^), and notice that the addition of the two Lie algebra generators in the two nearby points corresponds 
to a chord diagram in which a dashed line starts and ends in a given edge of the spin network. Another possibility 
is if the action is close to a vertex. In that case one can have the insertion of the two Lie algebra generators in two 
different links adjacent to the vertex. Concretely, we get. 




, K = (80) 



C{N)E ' Y ^K =0 (81) 



where we see that the contributions vanish in the first expression due to the contraction of the Cafcc with a symmetric 
expression in 6, c. In the second expression, given the action on an edge i, the contribution arising from connecting 
that edge with another edge j is equal and opposite to that arising from considering the action on the edge j and 
considering the connection with the edge i, and they therefore cancel each other as well. This is all only true because 
we have chosen a symmetric regularization, otherwise the expression would not be symmetric under the interchange 
of b and c. 

We should notice that this is a major departure from what we did in our previous paper p^ . There we did not 
choose the symmetric regularization, and the diffeomorphisms did not vanish on these invariants. At that point we 
interpreted the action of the constraint as associated with the change in the value of the framing invariants due to 
the addition of "writhes" by the diffcomorphism. It should be noted that there is some ambiguity in what one means 
by diffcomorphism in the case of framing dependent invariants, i.e., if one is talking of diffeomorphisms of "ribbons" 
or of loops with no width. Choosing the invariants to be diffcomorphism invariant is in line with adopting the point 
of view that one is considering diffeomorhpisms of "ribbons" . 

What we find here is that the motivation for the symmetric regularization that came from the connection representa- 
tion has a natural counterpart in loops: in the connection representation the regularization ensured that the constraint 
generate diffeomorphisms, and therefore that the Chern-Simons state be annihilated by the constraint. The same is 
true in the loop representation. With this regularization, the loop transforms of that state are diffcomorphism invari- 
ant. In other words, in this regularization, the diffcomorphism considers the framing dependent objects as invariants. 
All coefficients in the expansion of the expectation value of the Wilson loop, and in particular all the independent 
components of each coefficient (the Vassiliev invariants) we discussed before, are annihilated by the diffcomorphism 
constraint. That is, in the regularization chosen, both the framing independent and the framing dependent invariants 
are annihilated by the diffcomorphism constraint. It is still open to question if the use of framing-dependent objects is 
warranted in quantum gravity, where the whole formalism has been from the outset set up in terms of loops without 
any reference to a framing. 



V. CONCLUSIONS 



To conclude, we summarize on the results found up to now: a) We have a space of wavefunctions of spin networks 
that are loop differentiable (the Vassiliev invariants and all the coefficients of the expansion of the expectation value 
of the Wilson loops); b) The diffcomorphism constraint naturally constructed in terms of the loop derivative is well 
defined and annihilates these states. 

We have elucidated several aspects of the generalization of Vassiliev invariants to spin networks. To begin with, 
we see that in general the expectation value of a Wilson net cannot be factorized as a framing independent invariant 
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times a framing dependent prefactor as was the case for loops. This, however, is possible if the spin networks have 
trivalent or simpler intersections. It nevertheless possible to identify at each order in the perturbation expansion, and 
for a given spin network, invariant quantities that are ambient isotopic (framing independent) and we illustrate the 
procedure to find them up to second order in the perturbative expansion. We have given expressions for the invariants 
that are given by analytic integrals along the lines that form the spin network. This allows to operate with these 
wavefunctions in a direct manner with the operators we will discuss in the following paper. 

The fact that the operator we chose for the diffeomorphism constraint generates the geometrical action of diffeo- 
morphisms on Wilson nets, suggests that the algebra of constraints will reproduce the expected classical Poisson 
structure. At the moment however, we have only constructed states in the loop representation where the constraint 
vanishes and therefore the algebra is reproduced, albeit trivially. In the companion paper we will introduce "habitats" 
of non-diffcomorphism invariant states related to the ones we introduced here and we will explicitly show that the 
classical Poisson algebra is reproduced at the quantum level. 

Having a well defined action for the loop derivative opens the possibility of introducing quantum versions of the 
Hamiltonian constraint and to implement in a concrete well defined setting proposals for the Hamiltonian constraint 
that have been shown to have the correct algebra at a formal quantum level We will expand on this and other 
topics in the companion paper. 
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APPENDIX A: LOOP DERIVATIVE OF FRAMING INDEPENDENT VASSILIEV INVARIANTS 

In this appendix we will prove a property of a limiting behavior of the loop derivative of the framing independent 
Vassiliev invariants based on trivalent intersections. What we will show is that in the limit in which the path on which 
the loop derivative shrinks to zero, the loop derivative vanishes. The interest of this limit is that it is the one that 
appears in the loop derivatives that arise in the diffeomorphism and Hamiltonian constraints of quantum gravity in 
the loop representation that we will discuss in the forthcoming paper. 

To prove this statement, we start by defining the invariant J(s, k). 



E{s, k) = exp 



E{s,0) 



< Wis) > 



(1) 



J(s, k) 



(Al) 



that is, we extract the framing dependence in the prcfactor, something that is allowed for trivalent or lower intersec- 
tions. 

Using the recoupling identities for trivalent intersections and equations (wW and (uq), one gets, 



lim 

e-+0 



E{s,0) 



< W{s) > 



(1) 



(A2) 



and comparing this expression with the derivative of (Al), and bearing in mind that the loop derivative satisfies 
Leibnitz' rule, one immediately concludes that 



lime^ / d^yx.(y,o)A,fc(7r^)J(s,K) =0. 



(A3) 



Since the result is independent of k, it implies that all terms in the expansion of J(s, k) in powers of n is annihilated 
by the loop derivative as well. Moreover, since we could have carried this construction for any gauge group (the only 
elements we used was that recoupling identities for trivalent intersections are tantamount to extracting a prefactor, 
which is true for any group), we conclude that the topological invariants are annihilated by this limit of the loop 
derivative. 
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